arXiv:1506.01401v2 [gr-qc] 23 Feb 2016 


SPECTRAL ACTION MODELS OF GRAVITY ON PACKED SWISS 

CHEESE COSMOLOGY 
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Abstract. We present a model of (modified) gravity on spacetimes with fractal structure 
based on packing of spheres, which are (Euclidean) variants of the Packed Swiss Cheese 
Cosmology models. As the action functional for gravity we consider the spectral action of 
noncommutative geometry, and we compute its expansion on a space obtained as an Apol¬ 
lonian packing of 3-dimensional spheres inside a 4-dimensional ball. Using information from 
the zeta function of the Dirac operator of the spectral triple, we compute the leading terms 
in the asymptotic expansion of the spectral action. They consist of a zeta regularization of 
a divergent sum which involves the leading terms of the spectral actions of the individual 
spheres in the packing. This accounts for the contribution of the points 1 and 3 in the 
dimension spectrum (as in the case of a 3-sphere). There is an additional term coming from 
the residue at the additional point in the real dimension spectrum that corresponds to the 
packing constant, as well as a series of fluctuations coming from log-periodic oscillations, 
created by the points of the dimension spectrum that are off the real line. These terms 
detect the fractality of the residue set of the sphere packing. We show that the presence 
of fractality influences the shape of the slow-roll potential for inflation, obtained from the 
spectral action. We also discuss the effect of truncating the fractal structure at a certain 
scale related to the energy scale in the spectral action. 


1. Introduction 

1.1. Fractal structures in cosmology. The usual assumptions of isotropy and homogene¬ 
ity of spacetime would require that the matter distribution scales uniformly in space. Large 
scale violations of homogeneity were discussed, for instance, in [17], while the idea of a frac¬ 
tal distribution of matter, scaling with a fractal dimension D ^ 3, was suggested in |46j . 
More recently, a growing literature based on the analysis of redshift catalogs at the level 
of galaxies, clusters, and superclusters has collected considerable evidence for the presence 
of fractality and multifractality in cosmology. We refer the reader to the survey [53] for a 
detailed discussion, see also [28]. While there is no complete agreement on the resulting 
dimensionality, partly due to difficulties in the interpretation of redshift data in estimating 
co-moving distances, multifractal models in cosmology have been widely studied in recent 
years. Cosmological models exhibiting a fractal structure can be constructed, adapting the 
original “swiss cheese model” of m The resulting models are usually referred to as Packed 
Swiss Cheese Cosmology (PSCC), see [S] for a recent detailed survey. The main idea in the 
construction of swiss cheese models of cosmology is to have spacetimes that are locally inho- 
mogenous but appear globally isotropic and that satisfy everywhere the Einstein equation. 
In the original construction of PSCC models, in a region dehned by a standard Friedmann- 
Robertson-Walker (FRW) cosmology, several non-overlapping spheres are inscribed, inside 
which the mass is contracted to a smaller higher density region, hence creating inhomo¬ 
geneities. The solution inside the ball is patched to the external FRW solution along a 
surface with vanishing Weyl curvature tensor (which ensures isotropy is preserved). A swiss 
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cheese type model based on the Tolman metric was developed in [19], [50]. In Packed Swiss 
Cheese Cosmology models, a conhguration of such spheres is chosen so that they are tangent 
to each other and arranged into a higher dimensional version of the Apollonian packing of 
circles. In a variant of this model, see the discussion in §8 of [H], instead of compressing 
the matter inside each spherical region, at each stage of the construction process the matter 
is expanded to he along the spherical shell, so that one ends up with a model of gravity 
interacting with matter, supported on the resulting fractal. The point of view we follow 
in this paper is similar to the latter: we consider spacetimes that are products of a time 
direction and a fractal arrangement of 3-spheres (or of other spherical space forms). We 
develop a model of gravity on such Packed Swiss Cheese Cosmology (PSCC) models using 
the spectral action as an action functional for gravity. 

1.2. Spectral triples. In Noncommutative Geometry, the formalism of spectral triples ex¬ 
tends ordinary Riemannian (and spin) geometry to noncommutative spaces, [TB]. This ap¬ 
proach encodes the metric structure in the data of a triple ST = [A^TL^D) of an involutive 
algebra A (associative, but not necessarily commutative), with a (faithful) representation 
•n ■. A ^ B{'H) by bounded operators on a Hilbert space "H, and with the additional structure 
of a Dirac operator D, namely an unbounded, self-adjoint operator, densely defined on "H 
with the properties: 

(i) : (/ -|- is a compact operator 

(ii) : for all a E A, the commutators [D, 7r(a)] are densely defined and extend to bounded 
operators on H. 

The metric dimension of a spectral triple is defined as 

(1.1) dsT ■= inf{p > 0 I tr((J + < oo}. 

A spectral triple is said to be finitely summable if dsr < oo. 

The notion of dimension for a spectral triple is more elaborate than just the metric di¬ 
mension. Indeed, a more refined notion of dimension is given by the dimension spectrum, 
Tst C C. This is a set of complex numbers, defined as the set of poles of a family of zeta 
functions associated to the Dirac operator of the spectral triple. In the case where Ker D = 0, 
the zeta function of the Dirac operator is given by Cd{s) = Tr(|D|“®). Let S{T) = [|D|,T] 
and let B denote the algebra generated by the S"^{7i{a)) and 6^{[D,7i{a)]), for all a E A, 
and m G N. One considers additional zeta functions of the form C,D,a{s) = Tr(a|D|“^), for 
arbitrary a E A and CD,b{s) = Tr(6|D|“^), for arbitrary b E B. The dimension spectrum is 
the set of poles of the functions CD,a(s) and CD,b{s). It typically includes other points, in 
addition to the metric dimension, and may include real non-integer points as well as com¬ 
plex points off the real line. Spectral triples associated to fractals typically have non-integer 
and non-real points in their dimension spectrum. In the following, we will use the notation 
Tgrp := n M+ for the part of the dimension spectrum contained in the non-negative 
real line. Geometrically, the dimension spectrum represents the set of dimensions in which 
the space manifests itself, when viewed as a noncommutative space. Even in the case of 
an ordinary manifold, the dimension spectrum contains additional points, besides the usual 
topological dimension. The non-negative dimension spectrum T'gj,, in particular, describes 
the dimensions that contribute terms to the action functional for gravity, as we discuss more 
in detail in 1 11.31 below, while the points of the dimension spectrum that lie off the real line 
contribute fluctuations in the form of log oscillatory terms, as we will see in 1 13.11 We say 
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that the dimension spectrum is simple if the poles are simple poles. Spectral triples with 
simple dimension spectrum are sometimes referred to as “regular”. However, the terminol¬ 
ogy “regular spectral triple” is often used in the literature with a different meaning, related 
to “smoothness” properties (see for instance [IB]). Thus, in the following we will use the 
terminology “simple dimension spectrum” to avoid confusion. 

A compact spin Riemannian manifold M can be described by a spectral triple STm = 
IPm), by taking A = the algebra of smooth functions, "H = 

L^(M, §) the Hilbert space of square-integrable spinors, and D = IpM the Dirac operator, 
which is a self-adjoint square root of the (negative) Laplacian of the manifold. The metric 
dimension of STm agrees with the dimension of M, by Weyl’s law for the Dirac spectrum. 
One can also recover the geodesic distance on M from STm- for any two points x,y & M 

dgeo{x,y) = sup{|/(a:) - f{y)\ \ ||[D, 7r(/)]|| < 1}. 

A reconstruction theorem mi moreover shows that the manifold M itself can be recon¬ 
structed from the data of a commutative spectral triple that satishes a list of additional 
axioms describing properties of the geometry such as orientability, Poincare duality, etc. 
The non-negative dimension spectrum consists of non-negative integers less than or 

equal to dim(M) (see 1 12.91 for more details). 

1.3. The spectral action as a model for (modified) gravity. The formalism of spectral 
triples plays a crucial role in the construction of models of gravity coupled to matter based 
on noncommutative geometry. The main ideas underlying the construction of these models 
can be summarized as follows: 

• The spectral action is a natural action functional for gravity on any (commutative 
or noncommutative) space described by a finitely summable spectral triple. 

• On an ordinary manifold, the asymptotic expansion of the spectral action recovers 
the usual Einstein-Hilbert action of gravity, with additional modified gravity terms 
(Weyl conformal gravity, Gauss-Bonnet gravity). 

• In the case of an “almost commutative geometry” (locally a product M x F of an 
ordinary manifold M and a finite noncommutative space) the model of gravity on 
M X F given by the spectral action describes gravity coupled to matter on M, with 
the matter content (fermions and bosons) completely determined by the geometry of 
the finite noncommutative space F. 

We refer the reader to the detailed account of the construction of such models given in 
m and in Chapter 1 of [IH]. For a finitely summable spectral triple, the spectral action 
functional [S] is dehned as 

5(A) = Tr(/(B/A)) = Y1 Mult(A)/(A/A), 

AGSpec(D) 

where / is a non-negative even smooth approximation to a cutoff function and A is a positive 
real number. As A grows, more rescaled eigenvalues of the form A/A escape the cutoff of / 
and the expression grows. 

In the case of a hnitely summable spectral triple with dimension spectrum consisting of 
simple poles on the positive real line, the spectral action can be expanded asymptotically 
for large A, [Sj. The asymptotic expansion relies on the Mellin transform relation between 



4 


ADAM BALL AND MATILDE MARCOLLI 


the zeta function of the Dirac operator and the heat kernel. The asymptotic expansion of 
the spectral action is then of the form 

(1.2) Tr(f(D/A))~ ^ ffAl’+ f(0)(M0), 

where f{v) dv are the momenta of /, the summation is over the points of the 

non-negative dimension spectrum and the coefficients are residues of the zeta function, 

(1.3) = ^Res,=/3 Cd(s), 

representing the noncommutative integration in dimension fd. 

In the case of a 4-dimensional manifold M, one can write the asymptotic expansion in the 
form [9] 

Tr(/ {D/ A)) ~ 2A'^/4ao + 2A^/2a2 + foa^, 

where the fi are momenta of the cutoff function /, with /o = /(O) and fk = /q°° f{v) dv. 
Physically, the coefficients Oq, 02 and 04 correspond, respectively, to the cosmological term, 
the Einstein-Hilbert term, and the modihed gravity terms (Weyl curvature and Gauss- 
Bonnet) of the gravity action functional. In the case of an almost-commutative geometry, 
the asymptotic expansion of the spectral action delivers additional bosonic terms, including 
Yang-Mills terms for the gauge bosons, and kinetic and interaction terms for Higgs bosons, 
and (non-minimal) coupling of matter to gravity (with the Higgs conformally coupled to 
gravity). The fermionic terms in the action functional for gravity coupled to matter come 
from an additional term not included in the spectral action, which accounts for the kinetic 
terms of the fermions and the boson-fermion interaction terms, see HU, na. For the purpose 
of the present paper, we are only interested in the gravitational terms, though couplings to 
matter could also be included, by taking a product of the geometries we will be discussing 
with a finite noncommutative geometry. 

We will see in the next section that, in the case of the Packed Swiss Cheese Cosmology, 
the spectral action has new contributions that arise from an additional point in the dimen¬ 
sion spectrum that reflects the fractality of the model, as well as log-periodic oscillations 
contributed by the points of the dimension spectrum that are off the real line. 

1.4. Summary of results. The main new results in this papers are structured as follows. 

In 1 12.11 we obtain an estimate, in the form of an upper bound, on the exponent of con¬ 
vergence of the zeta function C£('S) of the length spectrum of an Apollonian packing V of 
3-sphere (Proposition [22]), we describe the spectral triple of V fDefinition l2.3p . and we com¬ 
pute the zeta function of the Dirac operator of the spectral triple (Proposition |2S|), 

in terms of the zeta function of the unit 3-sphere and the zeta function Cr:(^) of the length 
spectrum. We discuss the structure of the dimension spectrum fLemma 12.7p . 

In l|3l we use the results on the zeta function to obtain an expansion of the spectral action 
functional. In 1 13.1l we discuss how the heat kernel expansion, and consequently the expansion 
of the spectral action, is altered by the presence of complex points of the dimension spectrum 
off the real line. For the case of a fractal geometry with exact self-similarity realized by a 
single contraction ratio, we obtain an explicit form of the log-oscillatory terms coming from 
the non-real points of the dimension spectrum, in the form of a Fourier series that converges 
to a smooth function (Proposition [TT]). In 1 13.21 we discuss approximations by truncation of 
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the Fourier series of the oscillatory terms. We then identify a set of four analytic conditions 
on the zeta function C£('5) of the length spectrum of the Apollonian packing (Dehnition 
I3.3p . which ensure that the spectral action has an expansion where the oscillatory terms 
can be approximated by a series of contributions from length spectra (fractal strings) with 
exact self-similarity. The contribution from the real points of the dimension spectrum yields 
gravitational terms as in the case of a 3-dimensional geometry, with an additional term 
coming from the only real pole of C£('S) at its exponent of convergence (Proposition [33]) • We 
also compute the form of the expansion of the spectral action when taking a geometry that 
is a product of the Apollonian arrangement V of 3-spheres with a compactified time axis 
(Proposition 13.6p . 

In ^ we investigate the effect on the spectral action functional of a truncation of the 
fractal structure at a certain energy dependent scale. We obtain estimates on the size of the 
error term and its dependence on the energy A fPropositions 133 and 14.6]) . 

In §3] we construct another model of fractal space, which allows for the presence of “cosmic 
topology”. This is obtained by taking a Sierpihski fractal arrangement of spherical dodec- 
ahedra and then simultaneously closing up all of them via the action of the icosahedral 
group, obtaining a fractal arrangement of Poincare homology spheres (usually referred to as 
dodecahedral spaces in the cosmic topology literature). This is a simpler fractal than the 
Apollonian sphere packing, since it has exact self-similarity with a single contraction ratio 
(2 -|- 0)“^, where 0 is the golden ratio. In this case we can compute more explicitly the 
new terms that arise in the expansion of the spectral action, including the oscillatory terms 
(Propositions 15.21 and 15.31 and Corollary 15.4p . 

In §3] we compute the effect of the additional terms in the spectral action expansion on 
the shape of the slow-roll potential obtained by perturbing the Dirac operator by a scalar 
held (Propositions 16.11 and 16.2p . 


2. Spectral triples and zeta functions for Packed Swiss Cheese Cosmology 

2.1. Apollonian packings of D-dimensional spheres. Higher dimensional generaliza¬ 
tions of the Apollonian packings of circles in the plane, consisting of “packings” of {D — 1)- 
dimensional hyperspheres inside a D-dimensional space were variously studied, 

for instance in [26], [30], [36], [10], [13], [51]. We recall here some useful facts, following [30] . 

A Descartes conhguration in D dimensions consists of D -I- 2 mutually tangent (D — 1)- 
dimensional (hyper)spheres. We write for a sphere of radius a. The curvature c = 1/a 
is endowed with positive sign for the orientation of with an outward pointing normal 
vector and negative for the opposite orientation. The curvatures of the spheres in a Descartes 
conhguration satisfy the quadratic Soddy-Gosset relation 



This relation can be formulated in matrix terms as c^Qdc = 0, with c = (l/oi,..., 1 / 00 + 2 ) 
the vector of curvatures, and Qn the quadratic form determined by the matrix 

Qd — Id+2 D lo+2lo+25 
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where = (1,1,•••,1) and Id +2 is the identity matrix. The angmented curvature-center 
coordinates of a sphere with center x = {xi,... ,xd) in consist of a (H -|- 2 )-vector 

lIxlP — 1 1 1 

w = (-, -xi ,..., -xd), 

a a a a 

where the hrst coordinate describes the curvature of the sphere obtained from the given one 
by inversion in the unit sphere. The reason for the hrst coordinate is so that one can extend 
unambiguously the augmented curvature-center coordinates to include the special case of 
degenerate spheres with zero curvature (hyperplanes). Given a Descartes conhguration of 
spheres, one assigns to it a (D -|- 2 ) x (D -|- 2 ) matrix W whose j-th row is the vector of 
augmented curvature-center coordinates of the j-th sphere in the conhguration. The space 
Ain of all possible Descartes conhguration in D dimensions is then identihed with the space 
of all solutions W to the equation 

/ 0 -4 0 \ 

(2.2) >V*Qz)>V= -4 0 0 . 

yO 0 2IdJ 

The space of solutions is endowed with a left and a right action of the Lorentz group 

0{D + 1,1). 

The D-dimensional Apollonian group is the group generated by the (D -|- 2) x (D -|- 2) 
matrices Sj of the form 

c T 2 , 2 D , 

~ ^0+2 A _ ^^ 1^+2 ’ 

with Id +2 the identity, Cj the j-th standard coordinate vector, and li:i +2 the vector with all 
coordinates equal to one. 

It is shown in [30] that in dimension D > 4 the Apollonian group Qd is no longer a 
discrete subgroup of GL(D -|- 2,M) and its orbits on no longer correspond to sphere 
packings. However, the dual Apollonian group is a discrete subgroup of GL(D -|- 2,M), 
and the Apollonian packings of {D — l)-dimensional spheres we will be considering here are 
obtained, as in Theorem 4.3 of [30], as orbits of the dual Apollonian group on A4 d- The 
dual Apollonian group is generated by reflections Sj- of the form 

(2.3) Sj~ = Id+2 + 2 fD+2^\ — 4 Cje*, 

with Cj the j-th unit coordinate vector. The matrix Sj- implements inversion with respect 
to the j-th sphere. The Apollonian packing is obtained by iteratively adding new Descartes 
configurations of spheres obtained from an initial one by iteratively applying inversions with 
respect to some of the spheres. When D ^ 3 the only relations in the dual Apollonian group 
Qj^ are (Sj-)"^ = 1. Thus, the spheres added at the n-th iterative step of the construction of 
the Apollonian packing are in correspondence with all the possible reduced sequences 

Jfc 7 ^ jfc+i, V/c, 

acting on the point W & AAd that corresponds to the initial Descartes configuration. Glearly, 
there are {D + 2){D + 1)"“^ such sequences, hence the number of spheres in the n-th level 
of the iterative construction is 

n] = (D + 2)(D + 

In the following, we will focus on the case D = 4, of Apollonian packings of 3-spheres. 
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2 .2. Lengths, packing constant, and zeta function. We proceed as in [12], [T3| to 

associate a spectral triple to an Apollonian packing Vd of (H —l)-spheres in a dimensional 
space. As above, let 

(2.4) CDiVo) = {an,k, neN,l<k<{D + 2){D + I)"”'} 

be the list (with multiplicities) of the radii an,k of the {D + 2){D + 1)”“^ spheres that 
are added in the n-th stage of the iterative construction of the packing. 

The packing constant (or exponent of the packing), aoiVo) of a packing Vd of {D — 1)- 
spheres is dehned as the exponent of convergence of the series 


that is, 

(2.5) 


1D + 2 )( D + 1)’'-1 

E E 

neN k=l 


(D+2)(D+l)"-i 

(Td{Vd) = sup{s e M; : ^ ^ = cx>} 

nEN k=l 

(D+2)(Z)+l)"-i 

= inf{s e M; ; ^ ^ < oo}. 

nEN k=\ 


For s > aoiVn), one dehnes the zeta function (coi^) the sum of the series 

(D+ 2 )(D+l)"-i 

(2.6) a„(s) = E E “".S' 

n^N k=l 


The zeta functions the more general zeta functions of fractal strings considered 

in (3^, need not in general have analytic continuation to meromorphic function on the 
whole complex plane, but there are a screen S, namely a curve of the form S{t) + it, with 
S : M. ^ {—oo, CDiVo)], and a window W consisting of the region to the right of the screen 
curve S in the complex plane, where (cd (■5) tias analytic continuation. We refer the reader to 
[37] for a more detailed account of screens and windows for zeta functions of fractal strings. 


2.3. Packing constant and Hausdorff dimension. The residual set of an Apollonian 
circle packing consists of the complement of the union of all the open balls consisting of 
the interiors of the circles in the packing. It was shown in [1] that the packing constant a 2 , 
dehned as in fl2.5p is equal to the Hausdorff dimension of the residual set of the circle packing. 
In the higher dimensional setting the problem of estimating the Hausdorff dimension of the 
residual set of a packing of {D — l)-dimensional spheres is much more involved, but there 
are some general estimates, obtained in [33] and [32] . 

Consider the inhmum of the packing constants over all packings Vd, 

aD = inf (Td{Vd)- 
Vd 

Assuming all the spheres in the packing are contained in the unit ball , and denoting 
by the D-dimensional ball with dB^ = the residual set of the packing is given 

^n,k ’ jr o o 

by 


TI(Vd) = X 
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Let dimj 7 ( 7 ^(P£))) denote the Hausdorff dimension of the residual set and 

6d = inf dimH(7^(Pz))) 

'Pd 

the infimum over all packings of the Hausdorff dimensions. The upper entropy dimension 
h^{7l{VD)) of the residual set TZ{Vd) is defined as 

h+{n{VD)) = limsup-i^^^^^^^^^^, 

log e 

where for a set X, the number N^{X) counts the smallest number of sets of diameter less 
than 2e that cover X. The lower entropy dimension is dehned similarly, with a liminf instead 
of limsup. It is known that the entropy dimension provides an upper bound for the Hausdorff 
dimension. Then we have the following estimates ([3S] and |32jL 

Proposition 2.1. The radii an,k of a packing Vd satisfy ^ = 1 and ^ = oo, 

hence D — 1 < aniVn) < D- The infima satisy 6 d < o'd, and for individual packings 

dimHiniVo)) < h+{n{VD)) = cjd^Vd). 

The identity 'Yhn k ^nk ~ ^ follows from the packing property, namely the requirement that 
the residual set TZ{Vd) in the H-dimensional unit ball has zero H-dimensional volume. The 
value dimniVo) is not known exactly. Some estimates are obtained, with various methods, 
in [26], [13], [51]. We provide a simple rough estimate in 1 12.41 below, for the specific case of 
3-spheres. 


2.4. Dimension estimate. Let P = ^*4 be an Apollonian packing of 3-dimensional spheres 
^a„k' compute here a rough approximation to the packing constant (74 (P) of the Apol¬ 
lonian packing, defined as in fl2.5l) . 

Proposition 2.2. By replacing the collection of radii {an,k} in the n-th level of the Apol¬ 
lonian packing V of 3-spheres with a single value a„ = Nn/'jn, where jn/Nn is the average 
curvature in the n-th level, one obtains an approximate estimate of the packing constant, 

CT4,av{V) ~ 3.85193... 


Proof. As discussed above, the number of 3-spheres in the n-th level of the packing P is 
given by the number of reduced sequences in the generators of the group Qn, namely 

(2.7) Nn := : fixed n} = {D + 2){D + l )"-%=4 = 6 ■ 5"-'. 

Let 7 „ denote the sum of the curvatures of the spheres in the n-th level. 


( 2 . 8 ) 


In 


6-5"-i 


E 


1 


As shown in Theorem 2 of [ID], the generating function of the 7 ^ = jn(s) is 


(2.9) 


Gi)=4('u) 


(1 — x)(l — 4x)u 
1 — ~ 


where n = 70 is the sum of the curvatures of the P-l-2 = 6 spheres in a Descartes configuration 
that gives the level-zero seed of the recursive construction. We obtain from this an estimate 
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of the metric dimension by replacing the curvatures l/an,k with their averages over 1 ( 
We denote the resulting approximation to the dimension by a 4 ^avi'P)- This is given by 

log(6 ■ 5”-‘) 


levels. 


\ D-O'" " / 

series. Since the specihc value of u does not influence the large 

■h rr« I n /-\cr' TrAV» /> i - I n /~J tTT/~» r\ im 


We expand fl2.9p in a power series. Since tne speciuc value ol u does not mnu 
limit above, we look at the values for u = 1, and we obtain 


lx", 


n behavior in the limit 

Gd= 4 = y^7n(l) 

n=l 

^ (11 + \/TM)"(-64 + 9\/TM) + (11 - V^)"(64 + 9\/TM) 

~ 3" • 10 • 

This then gives a^^aviV) ~ 3.85193 ... as stated. 


□ 


2.5. A spectral triple on the Cayley graph of the dual Apollonian group. Let 7i 

denote the Cayley graph of the dual Apollonian group Qj^. Since for H 7 ^ 3 the group Q 


generated by the D + 2 reflections Sf of (12.3p . with the only relations of the form {SfY = 1, 
the Cayley graph Tb is an inhnite tree with all vertices of valence D + 2. We endow the tree 
Td with the structure of a finitely summable tree, in the sense of §7 of [12], by choosing a 
base vertex vq and endowing all the A„ = {D + 2){D + 1)"“^ edges at a distance of n steps 
from Vo with lengths i{en,k) = cin,k, equal to the radii of the spheres in the n-th level of the 
sphere packing. Then, as in Theorem 7.10 of [12] one obtains a hnitely summable spectral 
triple 


TT 


= {^To^'^To^'^Td) = ®eG£:(rD)(-TrB,77qe),T)£(e) + 

The involutive subalgebra Ar^ of the C*-algebra C{Td) is determined, as in [12], by the 
condition that / G A 77 , has [ 1 ) 77 ,, 7r(/)] densely dehned and bounded, where vr : C{Td) —t 
BfiHrD) is the representation by bounded operators on the Hilbert space of the triple. The 
pairs (^^^(e), 7l£(e)) are constructed as in the “interval spectral triple” of §3 of [12], with 
T-La = -h^([—a, a],/i) with the normalized Lebesgue measure + and Da with eigenvectors the 
basis elements 6 ^ = exp(i 7 rmx/Q!) with eigenvalue irm/a. The Dirac operator Dj-^ then has 
spectrum 

.7r(2m + 1) 


= { 


Spec(T>ri3) = { 

7r(2m + 1) 


e G E{Td), m G Z+} 


2a 


n.k 


2£(e) 

uGN, l<k<(D + 2)(D + 1)"~^, m G Z+j. 


The shift Ti/2i{e)I to the Dirac operator Dn^f.) is introduced in [T2] to avoid a kernel, so 
that the zeta function C,Drjy{^) = is well dehned. The zeta function of the Dirac 

operator of the spectral triple STj-^^ is given by 


Tr(|r>r„ 


OS+l 


7r° 


0C(s)Cn.(«), 


where ((s) is the Riemann zeta function, see §7.1 of [12]. The exponent of summability of 
the spectral triple (the metric dimension) is equal to the packing constant of fl2.5l) . 
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2.6. The spectral triple of a sphere packing. Suppose given an Apollonian packing 

Vd of {D — l)-dimensional spheres in We modify the construction above, by 

introducing the contribution of the {D — l)-spheres of the packing, through their 

respective spectral triples. We replace the data construction above, 

for an edge en,k of length i{en,k) = an,k, with new data of the form no-i), where 

°"n,k 

'HgD-i = S) is the Hilbert space of square integrable spinors on the {D — l)-sphere 

and V^D-i is the Dirac operator, with spectrum 

Spec(DgD-^i) = {A^,± = ± a-\ (-y- + ^) : ^ e Z+} 

and multiplicities 

Definition 2.3. The spectral triple of the Apollonian packing 

Pd = {Sfy : n e N, 1 < i < (D + 2)(D + 

is given by 

(2.10) i), 

i{e) £{e) 

where To is the Cayley graph ofQ^, as above, with edge lengths l{en,k) = ci,n,k, o-nd the data 
{Usp-uVsp- i) are defined as above. The involutive subalgebra A-p^ consists of f E C{Vd) 

^n,k ^n,k 

with [D-p^, 7r(/)] densely defined and bounded. 

The fact that this is indeed a spectral triple follows from the general results of [12] and 
[15] . In particular, the spectral action of the Swiss Cheese Cosmology model is obtained by 
considering the case of a packing of 3-dimensional spheres, 

(2.11) STp^c" := (AI-P4,'H-P4,'DP4). 

In order to compute the spectral action for the spectral triple of a packing of 3-spheres, 
we first recall some facts about the spectral action of a single 3-sphere. 

2.7. The spectral action on the 3-sphere. We start by recalling some very simple and 
well known facts about the round sphere and its spectral action functional. We will need 
these in the rest of this section as building blocks to construct the spectral triple and the 
spectral action for the Packed Swiss Cheese Cosmology. 

The Dirac operator on the 3-sphere with the round metric of unit radius has spectrum 
Spec(D53) = {n + ^} with spectral multiplicities Mult(n -|- |) = ?7,(?7, -|- 1), hence the spectral 
action takes the form 

(2.12) 553 (A) = Tr(/(D53/A)) = ^ n(n + l)/((n + ^)/A). 

Lemma 2.4. The zeta function of the Dirac operator is of the form 

(2.13) CD,,(s)=2C(s-2,|)-y,|), 
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where C{s,q) is the Hurwitz zeta function. The spectral triple STgs has simple dimension 
spectrum, with — {1)3}. The asymptotic expansion of the spectral action is corre¬ 

spondingly of the form 

(2.14) 553(A)~AV3-^A/i 

Proof. The result immediately follows by writing 

= J2^{k+ l)(k + 2)(k + |)-' ^J2^(k+ 

k>0 k>0 k>0 


The Hurwitz zeta function C(s, q) has a simple pole at s = 1 with residue one, hence Cvgsis) 
has simple poles at s = 1 and s = 3, respectively with residues ReSs=iCDg 3 (s) = —1/2 and 
ReSs= 3 Crig 3 (s) = 2. Then applying fll.2p . one obtains the spectral action expansion. In the 
constant term we have CDg 3 ( 0 ) = 2(/(—2,3/2) — (/(0,3/2)/2 = 0, since ({—2,3/2) = —1/4 
and C(0, 3/2) = —1. □ 

Corollary 2.5. In the case of a 3-sphere S'f with the round metric of radius a > 0, the zeta 
function is of the form 

(2.15) Cd,3 M = «‘(2C(s - 2, |) - lc(i>. |)). 

and the asymptotic expansion of the spectral action is given hy 

(2.16) 553 (A) ~(Aa)Vs-^(Aa)/i. 

Proof. The spectrum of the Dirac operator D 53 is a scaled copy ^{^-\-Z) of the spectrum of 
D 53 , and the multiplicities coincide. Thus, we have 


Tr(|Ds3|-) = f; 2n{n + 1) ( 

n=l ^ ' 


oo 

2a^ ^n(? 7 , + l)(n + 1 / 2 )“^ 

n=\ 


oo 

= 2a" ^(n + l/2)2(n + 1/2)"" 

n=\ 


~2 


E("+V 2 )- 

n=l 


oo s ^ 

= 2a* ^(n + 3/2)-'-=) - T + 3/2)- 

n=0 n=0 

When 5R(s) > 3 (the metric dimension of the 3-sphere), this simplihes to fl2.15p . 


□ 


A method for non-perturbative computations of the spectral action functional based on 
the Poisson summation formula was developed in |9], for sufficiently regular geometries for 
which the Dirac spectrum and the spectral multiplicities are explicitly known. In particular, 
the spectral action for the round sphere was computed in [9] using this method. The 
computation was generalized to spherical space forms, and to 3-dimensional tori and 
Bieberbach manifolds in HU, Ba, Ha. The computation of [9] for the 3-sphere can be 
summarized quickly as follows. Let / be a rapidly decaying even function. The eigenvalues 
of D 53 form an arithmetic progression, and there is a polynomial P(u) = ^ that 







12 


ADAM BALL AND MATILDE MARCOLLI 


interpolates the spectral multiplicities, Mult(A) = -P(A). Thus, one can write the spectral 
action as 

‘553(A) = + ^), 

where g{u) = — |)/(m/A) is also a rapidly decaying function. This is then the sum of 

values of a rapidly decaying function on points of a lattice, which can be evaluated using the 
Poisson summation formula 


n&, 


n -\—) = 

2 ' 




nGlj 


where 


g{x) = / = hu^- -)/(u/A)e-^"““du 

Jto -/to> 4 


is the Fourier transform 


1 


g{x) = A^ / v^f{v)e-^^^^^^dv - -A / f{v)e-^^^^^^dv, 

Jr 4 

after substituting u = Av. Let denote the Fourier transform of u^/(u), in the first term 
above. It is shown in |9] that the sum on the Fourier transformed side can be very accurately 
approximated by the term with u = 0 , yielding for any k 

Tr(/(D/A)) = [ v^f(v)dv-]A [ f{v)dv + 0{A-^). 

Jr ^ Jr 

In the case of the round 3-sphere of radius a, we have 

Tr(/(^)) = Tr(/(:^)), 

and the approximation formula above extends to S^, replacing A with Aa, so one obtains 


(2.17) 


Tr(/(T> 53 /A)) = (Aa)^ / f{v)dv - ^{Aa) I f{v)dv + 0{{Aa) ^), 


for arbitrary K E N, which agrees with the expression fl2.16p . with the error term as in [9]. 


2.8. Zeta function of a 3-sphere packing. We focus here on the case of a packing V = Vi 
of 3-spheres, where at the u-th iterative step in the construction one has 6 • 5"’“^ spheres, 
with radii Qn^k with k = 1,. .., 6 • 5”“^, starting with an initial Descartes conhguration of 
6 mutually tangent 3-spheres. As above, let £ = £4 = | u G N, A: G {1,. .., 6 • 5”“^}} 

be the length spectrum of the radii of all the 3-spheres in the packing. We consider the 
associated zeta function fl2.6p for D = 4, which we denote simply by Cc('S), 

(2.18) Co(s) - £ £ <.k- 

nG'M k=l 

Proposition 2.6. Let < 74 (P) be the packing constant of an Apollonian packing V of 3- 
dimensional spheres, as in fl2.5p . For s > a^{V), the zeta function of the Dirac operator V-p 
of the spectral triple ST psc of fl2.1ip is given by 

(2.19) H'(|I>nr*)= (2C(s-2,|)-lc(i>4))a(s), 
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where 'is the Hurwitz zeta function and Cc{s) is as in fl2.18p . 

Proof. Since 0 ^ Spec(i5 53 ), is invertible and so is then the Dirac operator V-p for the 
spectral triple STpsc- The metric dimension is then given by inf{/3 > 0 | Tr(|D-p|“^) < 00}, 
where the zeta function is given by 

00 6 - 5 "-i 

Tr(|D^r‘) = ^ E 

n=l k=l 

Each term in this sum can be computed as in (12.151) . We can then evaluate the zeta function 
of the spectral triple ST psci using the fact that the contribution of each sphere ^ is of 
the form TrdD^s ~ 2 , |) — |C('5, §)), and we obtain 

00 . X 00 

H'dDI-') = = I 2C(f -2,-) - 2C(f.2)) 

k=0 ' ^ ^ k=0 

for s > a, with a = max{3, 1 , cr4(P)} = (T4(V), where aDiVo), as in fl2.5p . is the packing 
constant of Vp, the exponent of convergence of the series Yh^nk- know from 1 12.31 that 
3 < (74 < 4, hence max{3,1, (J 4 ('P)} = a^ifP). □ 


2.9. Dimension spectrum. The dehnition of dimension spectrum we are using in this 
paper is the same as in [19]. It is slightly different from other versions in the literature, see 
[2T| and [33]. In particular, note that the dimension spectrum Em for an ordinary smooth 
manifold M of dimension n = dimM the dimension spectrum is given by the set Em = 
= {0; 1)2, •• • ,77,}, according to Example 13.8 of [20], or by Em = {ur G Z : m < n}, 
according to [T8|, p.22, and Proposition A.2 of [33]. The leading terms in the asymptotic 
expansion of the spectral action, which correspond to the gravitational terms in the action 
functional, arise from the points in E)[^ = Em H M+, which are the same in all cases, hence 
for our purposes the slight discrepancy between different versions of the notion of dimension 
spectrum adopted in the literature does not affect the results. 

In the following, we will focus on analyzing the poles in and off the real line of the 
zeta function Cx>p(s) = Trd'DpI”*). While these poles certainly contribute points to the 
dimension spectrum, there may, in principle, be additional poles coming from other zeta 
functions Cb,Vj,{s) = Tr(6|D-p|“*), for algebra elements b ^ B not equal to the identity. In 
the case of smooth manifolds, it is known (see for instance Proposition A.2 of [33]) that 
these zeta functions do not contribute additional poles. While there is no general result for 
arbitrary spectral triple, in the case of the spectral triple of a fractal geometry it is often 
suggested that the subalgebra of “smooth functions” should consist of functions that are 
supported on hnitely many levels of the fractal construction (for example, in the case of a 
Cantor set, that would mean locally constant functions). If the fractal is built out of pieces 
that are smooth manifolds (as in the case of a sphere packing) then one should also require 
that the functions are smooth on each smooth component. While this choice of smooth 
subalgebra does not necessarily have, in general, the same good analytic properties as the 
algebra of smooth functions on a smooth manifold, it is a natural choice in this setting. In 
this case, the fact that the additional zeta functions Cb,v-pi.s) do not contribute new poles 
can then be reduced to the known case of manifolds. When we discuss perturbations of the 
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Dirac operator by a scalar field, to obtain a slow-roll potential for inflationary models, we 
will assume that the scalar fields also live in this smooth subalgebra. 

The result of Proposition 12.61 then shows that the dimension spectrum of the spectral triple 
STpsc is given by the following set. 

Lemma 2.7. The dimension spectrum Tipsc of the spectral triple STp^^ consists of the 
union of the dimension spectrum of the 3-sphere, a single other real point a 4 {V), and a 
countable collection of points off the real line, lying in the window W where has analytic 
continuation. 

In general it is difficult to characterize precisely the position of the poles that are off the 
real line, except in the case of self-similar fractals. We will discuss how to obtain some 
control on the contributions of these points to the expansion of the spectral action in the 
following section. 

3. Spectral Action eor Packed Swiss Cheese Cosmology 

In this section we use the results of the previous section on the zeta function of the Dirac 
operator on the Packed Swiss Cheese Cosmology in order to study how the spectral action 
is affected by the presence of fractality. In particular, under some restrictive assumptions 
on the analytic properties of the zeta function C-c('5) of the Apollonian packing, and using 
the relation between the heat kernel and the zeta function and results on the asymptotic 
expansion of the heat kernel, we will obtain an expansion of the spectral action that contains 
the familiar gravitational terms of a 3-dimensional sphere, but also has additional terms 
determined by the residue of the zeta function at the packing constant, and a Fourier series 
of additional oscillatory terms coming from fluctuations produced by the presence of poles 
of the zeta function located off the real line. 

3.1. Zeta function, heat kernel, and spectral action on fractals. An asymptotic ex¬ 
pansion for the spectral action, in the sense of [S], is known to exist (see Theorem 1.145 of 
Hi) whenever there is a small-time asymptotic expansion for the heat kernel of the corre¬ 
sponding Dirac operator. In the case of an ordinary manifold, or an almost-commutative 
geometry, the heat kernel expansion is known by classical results on pseudo-differential oper¬ 
ators. For more general spaces, like fractal geometries, there are no analogous theorems that 
hold with the same level of generality, although several results on the heat kernel expansion 
on fractals are available, see for instance the detailed survey given in [22]. For some general 
results about Laplacians on fractals and heat kernels we also refer the reader to [Mj, |52j . 

Specifically in relation to the asymptotic expansion of the spectral action, cases where the 
zeta function has poles off the real line, which contribute log-oscillatory terms to the spectral 
action, were studied in [23] and [24] . 

The main new feature that arises in the case of fractal geometries is, as we have seen in 
the previous section, the presence of poles of the zeta function that are off the real line. In 
the case of the geometry of the Apollonian packings of 3-spheres we consider in this paper, 
those poles correspond to the poles off the real line of the zeta function C£('S) of the length 
spectrum C = C{V) of the packing. 

As discussed in §1-3 of [37], for general zeta functions of fractal strings C the distribution 
of the non-real poles can be very complicated. In the best possible case, which corresponds 
to fractals with a self-similar structure where the contraction ratios are all integer powers of a 
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fixed scale 0 < r < 1 (lattice case) the non-real poles lie, periodically spaced, on finitely many 
vertical lines. In cases with self-similar structure, but where the contraction ratios do not 
satisfy the lattice condition (non-lattice case), the poles off the real line have a quasi-periodic 
behavior and are approximated by a sequence of lattice strings. 

In the case of a length spectrum with exact self-similarity realized by a single contraction 
ratio r, the poles off the real line lie on the vertical line with 5R(s) = a, which is the Hausdorff 
dimension, and with periodic spacings of length i^gii/r) ’ s = a -|- with m G Z. 

We will discuss in l lS.ll an example of this kind, which is relevant to our cosmological models. 
In such cases with exact self-similarity, it is known (see §4 of [22]) that the contribution of 
the off-real poles to the heat-kernel asymptotic consists of a series of log-oscillatory terms. 
We have the following model case for this situation. 


Proposition 3.1. Let X be a fractal geometry with a Dirac operator Dx of the associated 
spectral triple with the following property: the eigenvalues of\Dx\ grow exponentially like 6"', 
for some b > 1, and the spectral multiplicities also grow exponentially like a” for some a > 1. 
Then the spectral action iSx(A) = Tr(/(Zlx/A)) has an expansion for large A of the form 

E 27rim 

where Sm = cr + and a = . The coefficients fs^ are given by integrals 

27r^Tn 

For sufficiently rapidly decaying test functions f{u), the Fourier series f^m coii- 

verges uniformly to a smooth function fa{d) of the circle variable 6 = modulo 27rZ. 


Proof. The zeta function has the form (£,^{s) = 'Yhn = (1 — with simple poles 

A logT’ with exponent of convergence a = The trace of the heat kernel 
has the “exponential form” 

(3.2) Tr(e-'^') = 

n 

for some constants a, &. Then it is known (see §4.2 of [22]) that one has a small-time 
asymptotics of the form 


(3,3) 


Tr(e 


-w 


log a 

2 t 2 log 6 


, log a Trim. , Trim , 

r(77TAT + — )e>‘P(-—log*) 


21og6;^ 21og6 log& 


log 6 


5^r(s^/2) 


2 log 6 

Indeed, through the Mellin transform relation between the heat kernel and the zeta function 


\D, 


r(i>/2) 




this corresponds to 

CdAs) = TldCxI-) = 5^ 


r(s„/2) 

r(s/2) • (s - Sm) ■ log 6 


-|- holomorphic 
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with poles at s = Sm = O' + with residue 1 / log b. To obtain then an expansion for 
the spectral action, one proceeds as in Theorem 1.145 of [19]. One considers a test function 
written as Laplace transform as k{u) = e~^'^h{x)dx, so that k{tDx ) ^ ^^^xh{x)dx. 

Using the expansion (13.dh one then has 

-,2 ^ r ( 5 m / 2 ) ^- s ^/2 


k(tDl) ~ Y. 


2 log 6 


X 


'"^‘^h{x)d. 


X. 


Since 3f?(sm) = cr > 0, we can write x as Mellin transform 

1 I 


X 


/2 _ 


r(Sm/2) Jq 


—ri) ±21 _i 
e n 2 ^ 


hence we obtain 


TUHtDl)) ~ 

m 

Then taking f{u) = k{v?) we obtain 




k{v)v ^ ^dv. 


k{v)v 


_1 


dv = 2 


f{u)u^^ ^du. 


We then set t = A ^ to obtain the form of the spectral action and the expansion 

E noo 

Aiog6 (/ f{u)u"^~^du) ReSs=smCDx{s), 

mGZ 


which gives (13. Ih . Using the relation between Mellin and Fourier transform, we can rewrite 
the coefficients 

fs^ = 7 ^ r —= f F(A)e-“dA = 27rF(-27rm), 

log b Jo M Jr 

where A = 1^ and F{X) = and F{^) = {27r)~^ J^F{X)e^^^dX is the Fourier 

transform. Provided the test function / is sufficiently rapidly decaying, the function F{X) is 
also a rapidly decaying function (at A —)■ +cxo because of the behavior of / and at A —)■ —oo 
because of the term Then the Fourier transform F{^) is also rapidly decaying, hence 
the Fourier series 22m ^ fsm = L^m converges to a smooth function fo-(d). □ 


More generally, in the case of exact self-similarity realized by a set of contraction ratios 
{ri,..., r^}, the zeta function C£('S) has a denominator of the form 1 —exponent 
of convergence is the self-similarity dimension given by the real number a satisfying the self¬ 
similarity equation ~ scaling factors satisfy the lattice condition, 

namely if the subgroup HJli 1® discrete, then (see Theorem 2.17 of [37]) the 

complex poles lie on finitely many vertical lines with 5R(s) < a, and are periodically spaced 
with period 27r/log(r“^), where r is the multiplicative generator of the scaling group, with 
Tj = for some integers nj. In this lattice case, assuming all the poles are simple and there 
are no cancellations from a numerator of Cc{s), one still obtains an asymptotic expansion 
of the form fl3.3p with one log oscillatory series for each of the finitely many vertical lines 
containing the complex poles of C£('S)- 

In the case with exact self-similarity realized by a set of contraction ratios {ri,... ,rm} 
that do not satisfy the lattice condition, it is no longer true that the complex poles he on 
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finitely many vertical lines. It is known (Theorem 2.17 of [37]) that in this case there are 
no other poles on the line 3?(s) = a except the real pole s = a, but there is a sequence 
of complex poles approaching the vertical line 3fJ(s) = a from the left. Moreover, all the 
complex poles are contained in a vertical strip (Tq < 3f?(s) < a, for some ao G M. Moreover, 
in this general non-lattice case, the complex poles can be approximated by the poles of an 
inhnite family of lattice cases, with increasingly large oscillation periods (see §3 of jSTj), 
which in turn correspond to an inhnite family of Fourier series of log-oscillatory terms. 

Remark 3.2. In the case of the D = 2 Apollonian circle packings, there are known results 
that characterize the presence of self-similarity, [15]: these packings correspond to quadratic 
irrationalities, via a continued fractions argument. However, analogous results for the higher 
dimensional Apollonian packings, characterizing the presence of exact self-similarity,s are 
not presently known. 


3.2. Approximations and expansion. In more general situations, even for nice cases of 
fractal geometries with exact self-similarity, we do not have such explicit control over the 
oscillatory terms as in the case of Proposition 13.11 where one has a single scale factor for 
self-similarity. In particular, in cases of self-similar geometries that do not satisfy the lattice 
conditions, the oscillatory terms can only be described via a sequence of approximations. 
Thus, we need to introduce some choices of approximations in the description of the log- 
oscillatory contributions to the spectral action coming from the poles of the zeta function 
that are off the real line. 


A hrst, very rough approximation, which we will occasionally use in the following, consists 
of replacing the smooth function fa{0) in the expansion iSx(A) ~ A®" fo-(6*(A)) of Proposition 
13.11 with its average value on the circle. This corresponds to selecting only the zero order 
Fourier coefficient 


1 


- 27 r 


U0)d9 = U = 


log (5) 


f{u)u'^ ^du = ReSs=aCDx{^) ' / ^du. 


This corresponds to only counting the contribution of the pole s = a on the real line and 
neglecting the contributions of the poles that lie off the real line. 

In a similar way, one can decide to approximate the function fo-(6*(A)) by truncating the 
Fourier series at a higher order. In the case of a fractal geometry with the non-lattice 
property, where there is an inhnite sequence of lattice approximations (§3 of [37]) to the 
off-real poles of the zeta function, these give rise to terms with increasingly long oscillation 
periods in the expansion of the spectral action. One can then choose to truncate the Fourier 
series at some hxed size M = m/\ogh, so that oscillatory series with longer oscillation 
periods get truncated earlier and contribute less to the approximation. 


Note that truncating the Fourier series in the spectral action expansion at some size 
M = m/ log 5 can also be seen as truncating the series of oscillatory terms in the heat kernel 
expansion 03.31) . The size of these terms is determined by the size of the value of the Gamma 
function F(^^^ -|- ^|). The Gamma function decays exponentially fast along the vertical 
line 3?(s) = hence these oscillatory terms in the heat kernel expansion become rapidly 

very small in comparison to the contribution of the m = 0 term. 


In all of these cases, when we introduce approximations to the oscillatory terms, the 
approximation we obtain for the spectral action is no longer really an “asymptotic expansion” 
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in the sense of |3T]. The usual meaning of asymptotic expansion implies that the function 
can be approximated around some value of the argument (or a limit value) up to arbitrary 
high order. For the purpose of building gravitational models, it will suffice to obtain an 
expansion of the spectral action up to order (including the oscillatory terms), and some 
sufficiently good approximation in cases where the oscillatory terms cannot be fully computed 
explicitly. For this reason, in the following we will simply use the terminology “expansion” of 
the spectral action, rather than insisting on the stronger properties of a genuine asymptotic 
expansion. 

3.3. Some analytic assumptions. As we pointed out in Remark 13.21 unlike the Apollonian 
circles case, in dimension D = 4 we do not have a characterization of the presence of exact 
self-similarity in the sphere packing. However, in order to obtain a reasonably behaved 
model, with respect to the properties of the zeta function and the spectral action functional, 
we restrict our attention to a subset of all the possible Apollonian packing, identihed by a 
set of requirements on the properties of the associated zeta function C/:('S)- 

Definition 3.3. A packing V of 3-dimensional spheres is analytic if it satishes the following 
properties: 

(1) The zeta function Cc{s) of the packing V has analytic continuation to a meromorphic 
function on a region of the complex plane that contains the non-negative real axis. 

(2) The analytic continuation Cr:('5) has only one pole on the non-negative real axis, 
located at s = a^iV). 

(3) The poles of C/:('S) are simple. 

(4) There is a family n G M, of self-similar fractal strings with the lattice prop¬ 
erty, and with increasingly large periods, such that the complex poles of C£('S) are 
approximated by the complex poles of Cz:„('5)- 

In terms of screens and windows, as in 133 , the hrst condition above consists of the 
property that the screen function S' : M —)■ (—oo,(T 4 (P)] satishes S'(O) < 0. 

In the last condition, the period of a self-similar fractal string £„ with the lattice property 
is the length 7r„ := with the property that all the poles of Cc^i^) off fh® I'oal line he 

on hnitely may vertical lines 3fJ(s) = aj with periodic spacing by ■ The approximation 

condition means that, for all e > 0, there exists an n G N an i? = R{e,n) > 0, such that, 
within a vertical region of size at most R, the complex poles of C/:('5) are within distance e 
of the poles of Ccn{s). For more details see §3.4.1 of [37] and see Figure 3.6 of [37] for an 
explicit example of such an approximation. 

3.4. Heuristics of analytic assumptions. At present, we do not have a characterization of 
the locus of packings satisfying the constraints listed in Dehnition 13.31 ffor example, in terms 
of a geometric locus in the conhguration space M.d oi Descartes conhgurations). We can, 
however, provide some heuristic explanation for the geometric meaning of the requirement 
that the zeta function Ccd{^) of length spectrum Co = {an,fc} of Apollonian packing 
Vd of {D — l)-dimensional spheres satishes these properties. 

Consider the possibility that a sphere packing has exact self-similarity. This would mean 
that there is a hnite set {ri,... ,rm} of scaling ratios, with the property that, for all n, k, 
the radii an,k ^ 1^+ of the packing belong to the subgroup ®+- Tf^i® subgroup 

will, in general, be dense in (non-lattice case). In such cases, for the zeta function of 
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fl2.19p . the factor C£('S) would have analytic continuation to a meromorphic function on all 
of C (see Theorem 2.4 of 1371 ), hence the hrst condition of Dehnition 13.31 would certainly 
be satished. Moreover, in such a case, the second condition would be satisfied by Theorem 
2.17 of [3Z]. The third condition would be satished, at least in the general case (again by 
Theorem 2.17 of [2Z])- The last condition is obvious in the lattice case, and is a consequence 
of the approximation result of §3 of [57] in the non-lattice case. In the non-lattice self-similar 
case, the self-similar strings £„ with the lattice conditions are constructed using Diophantine 
approximation (Lemma 3.16 and Theorem 3.18 of |37jL 

Thus, one should think of the conditions of Dehnition 13.31 as a generalization of the good 
conditions satished by the zeta function of a fractal with exact self-similarity. One can expect 
that they may be fulhlled by especially regular (especially symmetric) choices of Descartes 
conhguration, although we leave a more precise mathematical investigation of this question 
to future work. 

3.5. Spectral action expansion. For the rest of this section we make the assumptions 
that the packing V of 3-dimensional spheres we are considering satishes the conditions listed 
in Dehnition 13.31 

Under the assumptions of Dehnition 13.31 the result of Lemma 12.71 on the dimension spec¬ 
trum Y^sTpsc rehned to the following form. 

Proposition 3.4. For a packing V of 3-spheres satisfying the properties of Definition \3.3[ the 
non-negative dimension spectrum of the spectral triple of the Packed Swiss Cheese Cosmology 
consists of the points 

(3-4) = {1, 3, (T4(P)}, 

with the metric dimension of the spectral triple is dpsc = The spectral triple has 

simple dimension spectrum, which, in addition to the points of on the real line con¬ 

tains a countable family of points off the real line, contained in a horizontally bounded strip 
o'min < 3f^(’S) < (74(P), approximated by the non-real poles of a family Cn of self-similar 
fractal strings with the lattice property and with increasingly large oscillation periods. 

Proof. Under the assumptions that the zeta function Cc{s) has analytic continuation to a 
window including the positive real axis, we see that the zeta function (vpis) = TrdP^pl”^) 
of the Dirac operator of the spectral triple STpsc also has analytic continuation to a mero¬ 
morphic function in the same region. Moreover, the assumption that the points {1,3} are 
not poles of C£(s) ensures that Cvpi^) has simple poles at these points. It also has a simple 
pole at s = UifP) and at all the poles oh the real line by the third assumption. Thus, the 
spectral triple has simple dimension spectrum and the non-negative part of the dimension 
spectrum is given by fl3.4l) . The last property about the poles oh the real line follows from 
the last property of Dehnition 13.31 and Theorem 2.17 of [37] . □ 

Proposition 3.5. LetV be a packing for which the assumptions of Definition \S.3\ hold. Then 
the expansion of the spectral action for the spectral triple STpscifP) is of the form 

(3.5) Tr(/(B„/A)) ~ A’ (c{3) A-A1 (c)!) /. + A’ (((a-2, 1) - 5 )) K, f,+SrW, 

where a = CilfP) the packing constant, 71^ = R.eSs=o-C/:('5) the residue of the zeta function 
of the fractal string C = jC{V), and fg = f{v)dv, the momenta of the test function, 

and Sp'^{A) is an oscillatory term involving the contributions of the points of the dimension 
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spectrum that are off the real line. For R > 0, let be the truncation of the 

oscillatory terms that only counts the contribution of the off-real poles with |Q'(s)| < R. 
Then the oscillatory term can be approximated by a sequence 

Nn 

(3.6) SrW<n ~ 

3=0 


where n —)■ cxo as R oo, and where anj = for 

_ 2'Km 

OTn,! “ 1 “ “ 1 “ 7 ) S j=0,...,Nn,m(3^ 

hgbn 

the set of non-real poles of the zeta functions Ccni^); CTmm < < cr and periods 

2ti/ \oghn —)■ oo and n ^ oo. The (6*„(A)) are smooth functions of the circle variable 
On = with Fourier expansion fa„j( 6 *n) = Em with 

3 1 3 

fsn,j,m ~ (C(^nj,m ~ 2, —) — —C{Sn,j,mi (s) / f (u) U du. 


Proof. Under the assumptions listed in Definition [33] on the zeta function Cn(s)) the residues 
at the points s = 1 and s = 3 of the dimension spectrum are given, respectively, by 

—13 1 

Res,=iC 7 ?p(s) = —Res,=iC(s, -) • Cz:(s) = -- C/:(l) 


Res,=3Ci?p(s) = 2Res,=3C(s “ 2, ^) ■ Cn(s) = 2C£(3). 

Thus, the terms in the expansion of the spectral action 

Tr(/CD„/A)) ~ LA'>/|Cp|-'’, 

h^^STpsc 

w\th.-^\V'p\~h the residues, as in (II.3p . are given by 

poo 

Tr(/(Dp/A))~ A^Cn(3) / v‘^f{v)dv 

Jo 

- Ajcai) r/(L<i<' 

3 13 

+ A^(C(ct- 2,-) --C(a,-))Res,=^C£(s) J v'^-^f{v)dv, 

where a = <74 (P). The approximate form of the oscillatory term is derived from the last 
property of Definition 13.31 and from the form of the oscillatory terms of Proposition 13.11 In 
the case of a self-similar string £„ with the lattice property, the poles of Ccn{s) off the real 
line consist of a finite union of sequences of the form Sn,j,m = (Tnj + i{o!n,j + with 

the same period and with an,o = 0 and anj 7 ^ 0 for j > 0 , and with an,o = o'n the 

self-similarity dimension, as shown in Theorem 2.17 of (3^. The terms in the expansion of 
the spectral action that correspond to these poles are then approximated, for large n, by a 
finite sum of terms as in Proposition 13.11 □ 
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For the purpose of this paper we will not give a more detailed analysis of the convergence of 
the approximation by the sequence A more precise analytic discussion 

of the nature of the approximation in (Id.bh will require a more detailed understanding of 
self-similar structures in higher-dimensional Apollonian sphere packings than is presently 
available, and will need to be addressed elsewhere. In terms of the expansion of the spectral 
action we are going to use in explicit gravitational models, we will truncate the series of 
oscillatory terms as discussed in 1 13.21 


3.6. Zeta regularization. The expression 03.51) for the spectral action of the sphere packing 
should be regarded as a “zeta regularized” form of the divergent series 

CX) 

5p(A) = 5^5s|.(A) 

k=0 

that adds the contributions coming from the spectral actions of the individual spheres in the 
packing. Indeed, since the spectral action of an individual sphere is of the form 

(3.7) ‘^53^ = A^a^ fc/s - ^Aa„,fc/i + 0((Aa„,fc)"^), 

and both points 1 and 3 are smaller than the exponent of convergence of the series 

(3.8) AVs - ^A/i ^ 

n,k n^k 

is divergent and requires a suitable regularization. The spectral action 03.5p can be inter¬ 
preted as such a regularization. Notice also that the error term 0{{Aan,k)~^) is very small 
for a fixed radius an,k and for sufficiently large A, but when the radii an,k vary over the set 
jC{V) of lengths of the packing V it becomes large for any given A, so that 03.7p cannot be 
extended directly to the whole packing. The term 

(3.9) A3/3C/:(3)-^A/iC£(1) 

in 03.5p is just a classical form of zeta regularization of the series 03.8p . with the divergent 
J2^n,k replaced by C£(3) and the divergent '^an,k replaced by Cz:(l). Th® additional term 
in 03.5p . which depends on the residue of Cz;(s) at s = (T^iV) detects the presence of a fractal 
structure in the geometry. We discuss these issues further in iQ] below. 


3.7. Packed Swiss Cheese Spacetime spectral action. In the previous subsection we 
computed the expansion of the spectral action for an Apollonian packing of 3-spheres, under 
some assumptions on the behavior of the associated zeta function. Here we consider an 
associated (Euclidean) spacetime model. This generalizes to the case of a packing of spheres 
the simpler case of a single sphere S^, where the associated spacetime is just M x S^, with 
the Euclidean time line M compactihed to a circle of size f3. We generalize the form of 
the spectral action of x Ff, by replacing the 3-sphere with a packing V of 3-spheres 
^, and using the results in the previous section. 
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Proposition 3.6. LetV he a packing of 3-spheres satisfying the three conditions of Definition 
\3.3[ Consider the product geometry Sjj x V of V with a circle of size (3. Then the spectral 
action has expansion with leading terms of the form 

(3.10) 

5s.,p(A)~ 2/l('A‘‘Ca3)l)3-A2lCc(l)lli + A''+' U(<y - 

where a = cr^ifP) is the packing constant fl2.5l) . T^o- is the residue of C,c{s) at s = a, and 


(3.11) 


f)3 := TT / h{p )p dp, f)i := 27r / h{p )pdp. 


(3.12) 


i)a = 2 h{p )p^dp. 


The oscillatory contributions from poles off the real line are approximated by a seguence 
(3.13) Ss^^MA)Tr ~ I] A‘'-..+‘t,,.,^(9„(A)), 

j=0 


with n ^ oo when R —)■ oo, where an,j = with Sn,j,m the non-real poles of 

and .{9n{A)) smooth functions with Fourier coefficients 

3 1 3 

with i)sr,jm defined as in fl3.12p . 

Proof. As observed in Lemma 2 of [9], the spectral action for Si x S'f, with the Dirac operator 


^Six53 - 


0 

D53 (g) 1 - i (8 Dgi 


D53 ® 1 + i ® Dsi 
0 


is of the form 

(3.14) 


TKAlflj.sj/A)) ~ 2/3ATr(K(D|./A)), 


for a test function h{x), and with the test function k on the right-hand-side satisfying 
A((r^) = /jg h(x^ -|- y‘^)dy. It then follows that the expansion of the spectral action on Sjj x S'f 
is given by (see Theorem 3 of [9]) 

Tr(A(Bjj,s3/A)) ~ 2/3 tA^o* fe - A"of, A 

with the notation of fl3.1ip . We now consider a similar situation, with the product geometry 
Sjj X Si replaced by Sj^ x V, where P is a packing of 3-spheres satisfying the conditions of 
Dehnition 13.31 The Dirac operator of the product geometry is again of the form 


D 


sUv 


0 D-p (g) 1 -j- ^ (g ^s\ 

Vp 0 1 — i g Dgi 0 
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where V-p is the Dirac operator of the spectral triple STpsc described in Dehnition 12.31 The 
same argument as in Lemma 2 of [9] shows that, as in fl3.14l) 

lV(ft(D|,,„/A)) ~ 2/3AlV(^(PpA)), 

with the test functions h and k as above. Using the result of Proposition 13.51 we then obtain, 
as above, the expression (13.101) . with [)„ given by 

p poo poo 

l)o-= / h{x‘^ + y‘^) dx dy = / h{p^)p'^dp / cos{9) d9 = 2 / h{p^)p''dp. 

J'Rj^x'R Jo J—nl2 Jo 

The structure of the oscillatory terms is obtained as in the previous Proposition. □ 


In cosmological models based on the spectral action (see mi, 021), the parameter fd is 
an artifact introduced by the choice of a compactification of the Euclidean time coordinate 
along a circle of size (d. As discussed in §3.1 of [12], the parameter fd can be interpreted as 
an inverse temperature and related to the temperature of the cosmological horizon. 


4. Fractality scale truncation 


A realistic model of fractal structures in cosmology will necessarily involve a choice of 
scale at which fractality is cut off: while the universe may involve a fractal structure at the 
scale of galaxy superclusters and clusters, it does not appear fractal at our scales, hence 
the self-similarity property is expected to break down at some level. In a gravity model 
based on the spectral action, which already naturally involves a dependence on an energy 
scale A, it is natural to assume that the scale at which fractality breaks down will be in 
some way dependent on A. In the construction of the spectral triple of the PSCC model, 
discussed in 1 12.61 above, we obtained a spectral action functional as a suitable kind of “zeta 
regularization” of the divergent series 


OO Nn 




n=0 k=l 


where the sum is over all the 3-spheres in the packing V, with an,k their radii, and with 
Nn = 6 • 5”“^, the number of spheres in the n-th level of the packing construction. Indeed, as 
we have seen in the previous section, the spectral action iS 7 [>(A) involves a zeta regularization 
of the above series, given by 

( \ reg 

0.1,k - On,k j = A^/s C£(3) - -A/i C£(l) 

n,k n,k / 

and an additional term 

A‘’(CU-2,|)-icu|))K,/„ 

involving the residue Ha = ReSs=o-CT(s) at a = (T 4 (P), which describes the fractality of the 
Apollonian packing. 
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4.1. Sphere counting function. In a model where fractality is truncated at a certain scale, 
one only considers the sphere packing V only up to a certain size. This requires estimating 
the number 


(4,1) 




of spheres in the given packing whose radii are of size at least a. In the case of Apollonian 
packings of circles and of 2-spheres it is known, by a result of [5] , that the Hausdorff dimension 
of the residual set of the packing is equal to 


dim^,(7^(P)) 


lim 

a—>-0 


\ogMa{V) 

logo 


so that, for a —)■ 0, one has A/'„(P) ~ a dim_ff{7ecp))+o(i)^ was proved in [35] that, 

in fact, one has the stronger result AfaiV) ~ A general heuristic argument 

for the existence of a power law governing the behavior of the sphere counting function for 
sphere packings in arbitrary dimension is given in [1]. Let S(V) denote the exponent of the 
power law, so that, for a —)■ 0 


(4.2) 


A4CP) ~ 


In fact, the result of [5] shows, in the case of an Apollonian packing of circles, that ^{V) is 
equal to the packing constant <72 (P), which combined with the result of |1] then gives the 
identification with the Hausdorff dimension. The general argument of §2 of [5] is independent 
of the dimension an it shows that, in general, one has the estimate 


a-i-O log a 


Thus, if the sequence has a limit, then the limit has to be the packing constant ud{'P)- 


4.2. Spectral triple with truncation of fractality scale. Thus, in a cosmological model 
where fractality is truncated at a certain size a, one would consider a spectral triple of the 
form 




n^k : 


n,k : CLri,k^^ 


where Pq C P is the part of the packing that includes only those spheres ^ with an,k > 
and C C(Va) that satishes the bounded commutator condition with the Dirac operator. 
Correspondingly, in this case, which involves only hnitely many spheres, the spectral action 
would be of the form 


(4.4) 




n,k : 


Lemma 4.1. Let V be a packing of 3-spheres satisfying the properties of Definition \3.3l and 
with the property that the function F{a) = — log(A/'o,(P))/log(a) has a limit for « —)■ 0. 
Then the spectral action S'p^{K) diverges at least like when a —)■ 0. 

Proof. Each sphere contributes to the spectral action a term of the form 
(4-5) A^s an,k + 0{{Aan,k)~^)- 

n,k n,k 
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Using the power law fl4.2p for a —)■ 0 we estimate 

5^ alj, > and o„,i > a'-'OT+'-d). 

n,k n,k 

By (14.31) and the hypothesis on the existence of the limit, we know that S{V) = a 4 {V), which 
we know satishes (T 4 {V) > 3, so that the exponents above are negative. Thus, for a hxed A 
and for a —)■ 0, the sum 

n,k:an^P;.>a \ n,k n,k / 

diverges at least like the dominant term □ 


4.3. Truncation estimates on the spectral action. One then expects that there will be 
a good approximation to the spectral action iS;p^(A) of the PSCC obtained by truncating 
fractality at a certain A-dependent scale. We discuss here possible choices of a function 
a = a (A) that retain the property of having a good control on the error term of the spectral 
action iS-p^^^j(A). We hrst recall how one obtains explicit estimates for the error term in 
fl2.17p of the spectral action on a 3-sphere, for a particular class of test functions, following 
the argument of Corollary 4 of [9]. 

Lemma 4.2. In the case of a cutoff function of the form f{x) = P{Tix^)e~'"^^, where P is a 
polynomial of degree d, the error term e(Aa) in the spectral action 

‘553(A) = (Aa)^ [ v^f{v)dv -^{Aa) [ f{v)dv + e{Aa) 

Jr 4 Jr 

satisfies the estimate 

(4.6) |e(Aa)| < {Aaf{5 + 7d + 

whenever Aa > -\/((i -I- 1)(1 -I- log((i -|- 1)) and Aa > 1. The coefficient Cq to he the sum of 
the absolute values of the coefficients of the polynomial Q, with /(x) = Q{7ix‘^)e~'^^^. 


Proof. As in Corollary 4 of [9], we have f^‘^\x) = {x^Zi^nx^) -|- Z2(7rx^))e“™^, where Zi = 
—Q + 2Q' — Q", and Z 2 = ^(Q — Q'), while generally for x > 3k{l -|- \ogk). 

For n ^ 0 one then has 

f{nAa) = Q(7r(nAa)2)e-5(”^“)' • 


because, by hypothesis, nijiAafi > 3d{l -|- logd). Since the decay is more rapid than simply 
exponential, we can see that 


25:' 

n=l 


(nAa)' 


< 2e"^ 2 ^ g-f n(AaF < 2e-f + 2- 


g-f (Aa)2 


n=4 


-1-^62 


(Aa)2 


■5(Aa)^ /oq 

< + 2— - < 2.023e-f 

-1 + 4.8 “ 

where we used the assumption that Aa > 1. Thus, we have 

^/(nAa) < 2.023C'Qe-5(^“)'. 

n^O 
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Similarly, for we get 




V P\nAa) < 2.023(2 + 3d + -d'^)CQt 

TT 

n^O 

By looking at the series for the spectral action after applying the Poisson summation formula, 
we see that the above terms contribute to the error as 

|e(Aa)| < (Aa)32.023(2 + 3h + 

TT 4 

which can then be estimated from above as in fid.op . □ 

We can then adapt the error estimates of Lemma [4.21 to the PSCC model, by performing 
a truncation on the fractality scale, dependent on the energy scale A. 

Remark 4.3. Under the assumption that an^k > the error term in fl4.5p is at most 
0((Aa)“^). Thus, it is natural to consider a model where the cutoff of fractality should 
happen at a scale a related to A by the property that a (A) ■ A grows like a positive power 
of A, so that one maintains a good control on the error term for large A. 

Proposition 4.4. Consider a truncated packing Va of 3-spheres with an,k > where 
a = a(A) = for some 0 < 7 < 1. Then the spectral action, computed using a test 

function of the form f{x) = P{'Kx^)e~'"^ with P a polynomial of degree d, satisfies 

1 


5- 


'^A-1+7 


(A) = 


E 


*'n,k 


AVa- 


'^n,k 1 A/i + e(A) 


^a„,fc>A-i+7 

where the error term satisfies 




>A-i+7 


|e(A)| < A=>+”<‘-’)<„(5 + 7d + 

where Omax = iLiax{an,A;} is the largest radius in the packing V, and a is the packing constant. 
Proof. The estimate follows immediately from the previous Lemma, since we have 
|e(A)|< ^ A3a3;,(5 + 7d + d2)C'Qe-^A«n,U^ 


ar...*:>A-i+7 


< VA-.„(P)A=a5,„(5 + 7d + d=)Coe-i* ’ = A3+‘'<‘-T)a^„(5 + 7d + d=)Cge-i 




□ 


We consider a further possible way of truncating the spectral action, for spheres with 
smaller radii and the behavior of the error term. We start with the following error term 
estimate, for a single sphere. 

Lemma 4.5. For a test function of the form f{x) = P{Tix‘^)e~'^^^, for some polynomial P 
of degree d, and for MAa the truncation scale, consider the truncated sum for the spectral 
action, 

Ss».m(A)= 5^ Mult(A)/(A)= ^ (a2_1)/(,L), 

X<MAa X<MAa 


'Aa' 


Let N = inf{?7 , G | + Z | tt, > max{MAa, |}}. Then, assuming that 

Aa > + 1)(1 + log((i + 1)) and Aa > 1, 
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< -(Aa)"(l + -(Aa)")Cpe-i' 
l\a / TT TT 


the error term satisfies 

KAa)|=2^(A= 

\>N 

Proof. Since | is a strictly positive quantity when evaluated at the integers, its absolute 
value is always less than vf. Furthermore, all terms in the series are positive, so that we 
can work with vf instead oi vf — \ and our conclusions will remain valid. As in Lemma 
1121 we use the fact that x^e 2 < 1, for all x > 3/c(l + log/c). Then, as long as x > 
+ 1)(1 + log(d + 1)), we know that x‘^f{x) < . 

Let our point of truncation be MAa, so that our sum is over u < MAa. Let N be as in 
the statement. Then we have 

2 


X>N 


.(Aa)|^ 25 :(V-L/(T)< 25 :(AaL(£) /(^ 


Aa 


\>N 


A 


Aa 


X 


Aa 


< 2 


\>N 


{Aafi—Cpe 2 (Aa) < —(Aa)^C 


00 

yp e ^ Oa) 


{AafiCpe 2 (Aa) 


TT 1 
' 2 (Aa)2 


< -iAayCpe-2 


^ 

v=0 
2/^ 


TT 


TT 1 

' 2 (Aa)2 


_ TT 1 

Now we check that (1 — e 2 (aS^)-i < 1 + |(Aa)^. Dehne g{x) = . In the case x > 1, 

we have 


p-iA _ 


n=0 


(-l)"x 


n\ 


-n 1 

= 1 — x~^ + 


E 


X 


l-2n 


X 


—2n 


„^V(2n-l)! (2n)! 


= 1— X ^ + -X ^ — ^ (x — 


n=2 


X 


—2n 


2nJ {2n-l)\ 


< 1 — X 


-1 


-x 


-2 


since every term in the third sum is positive. In the case 0 < x < 1, substituting u = 1/x, 
with M > 1, we have 

1 11 + 2e-“ 1 + 2e-“ 


9{x) = 


< 1 + 2e" 


1 - e-“ 1 - e-“ 1 + 2e-“ 1 + (e-“ - 2e-2“) “ 

Compare this with 1 + 1/m: the term 2e““ decreases faster than 1 / m , and 2e“^ < 1 , so for 
all M > 1 , we have 1 + 2e““ < 1 / m . Thus, for all 0 < x < 1 we have g{x) < 1 + x. Together, 
these two cases give the desired inequality 


(1-e < 1 + A(Aa)2, 

TT 

which then gives the stated estimate for the error term |e(Aa)|. 

For a given scale M, we consider a summation as above, of the form 


□ 




E 


2 


2 (A^ - 7 ) / 


A 


Aai 


We are interested now in the case of smaller spheres, with an,k < a, where a = a (A) was a 
previously chosen cutoff. Thus, we relate the scale M to A accordingly, by assuming that 
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M = M(A) = a(A)A = A^. This, in turn, gives the lower bound ^ = |A which means 
that we are considering spheres with an,k > Q!(A) = 


Proposition 4.6. For a(A) = A and d(A) = |A d+7)^ Qjid sufficiently large A, we 


have 


rAi+T'a„,fc]-l 


Sf.<»,(A)= E E 2(V-7)/ 


f A-T'<Aa„,fc<AT' A=f 


A 

Aai 




where 


h(A)| < (^0 ^CpA2^+(^+'^)"e-5^'E 


Proof. For large A we can estimate the number A/'o(a)(P) of spheres with On^k ^ d(A), with 
d(A)“°' = (1)'^ where a = cr^iV) is the packing constant. For each sphere in the 

summation, we can apply the error term estimate of the previous lemma, with Aan,k E A'^ 
and we obtain 

|e(A)| < ^ 2^2 ^2,^ 

\2 J TT TT 

□ 


We can view the estimate of the error term of this sum as a way to estimate the effect of 
changing the cutoff scale from a(A) = A“^+'>' to Q!(A) = 


5. Related fractal models 

5.1. Fractal dodecahedra and cosmic topology. In [7], [TT], [12], cosmological models 
based on the spectral action functional of gravity are constructed for (compactihed. Eu¬ 
clidean) spacetimes of the form x E where Y is either a spherical space form or a flat 
Bieberbach manifold, and it is shown that the spectral action detects the different cosmic 
topologies through the shape of an associated slow-roll inflation potential. In particular, it 
is shown in [7], j4T], [54] that the spectral action for a spherical space form Y = S^/T is 
given by 

Sk(A) = Ess.(A), 

independently of the spin structure (even though the Dirac spectrum itself changes for dif¬ 
ferent spin structures). Of particular interest for cosmic topology is the case where Y is 
the Poincare homology 3-sphere (dodecahedral space), which is regarded as one of the most 
promising candidates for a non-simply connected cosmic topology, [H], [3U] . 

We consider here a different possible model with fractal structure, where the building 
blocks are spherical dodecahedra, folded up to form Poincare homology spheres, arranged 
in a fractal conhguration that generalizes the Sierpihski fractal to dodecahedral geometry. 
Other similar constructions can be done using other solids, and adapted to the other spherical 
space form candidates for cosmic topologies. These fractals are much simpler in structure 
than the Apollonian sphere packing described in the previous section, as the successive levels 
of the construction are all obtained by repeatedly applying the same uniform contraction 
factor. This makes the computation of the Hausdorff dimension immediate, as well as its 
identihcation with the exponent of convergence of the relevant zeta function. Moreover, it 
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is also immediately clear that an analytic continuation exists to a meromorphic function on 
the entire complex plane, hence all the properties can be checked more easily. 

More precisely, let Ya = S^/Xi 2 o be the quotient of a round 3-sphere of radius a by the 
isometric action of the icosahedral group Xi 2 o- The choice of a fundamental domain, given by 
a spherical dodecahedron, and the action of the group X 120 determine a tiling of consisting 
of 120 dodecahedra. The quotient 3-manifold Ya is a Poincare homology sphere, of volume 
Vol{Ya) = Vol{S!/Iuo) = gaX 

Consider now the following well known construction of a Sierpihski type fractal based on 
the dodecahedron. Starting with an initial (solid) regular dodecahedron, one replaces it with 
20 new regular dodecahedra, contained inside the volume bounded by the initial one, each 
placed in the corner corresponding to one of the 20 vertices of the original dodecahedron. 
Each of the new dodecahedra is a copy of the original one scaled by a factor of (2 -|- 0)“^ 
where 0 is the golden ratio. One keeps iterating this procedure on each of the dodecahedra 
in the new conhguration. Let VY,n be the union of dodecahedra obtained at the n-th step 
of the construction, where we simultaneously perform the identihcation of faces, in each 
dodecahedron, according to the action of X 120 , so that each is folded up into a Poincare 
homology 3-sphere. Let Vy denote the resulting limit, which in the set theoretic sense is 
given by the intersection Vy = ^n>iVY,n- The following fact then follows immediately. 


Lemma 5.1. The Hausdorff dimension of the resulting set is 
(5.1) 


dim^Xy = . = 2.3296... 


log(2 + (f) 

This is equal to the exponent of convergence of the series 

a(7V)W = 20“(2+ 


n>0 


which has analytic continuation to a meromorphic function on C, 

Cn(Py)(s) - 1 _ 20(2 + 0)-* 

with simple poles at the points 


(5.2) 




log(20) 27Tim 

log(2 + 0) ^ log(2 + 0) ’ 


with m G Z, 


all with the same residue 
(5.3) 


^GSs=s^Cc{Vy){s) 


1 

log(2 + 0) ■ 


We then construct a spectral triple for the conhguration Vy using the same procedure as 
in [12]. Let {C{Ya)^'HYa^ Dy^) be the spectral triple of Ya with "Hy^ = LfiYa.'B) the square 
integrable spinors and Dy^ the Dirac operator. See [41], [54] for a more detailed discussion 
of these data. 


Proposition 5.2. The spectral triple 
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with A-Py C C(Vy) satisfying the bounded commutator condition, and with a„ = a(2 + 0) ", 
has zeta function 

(5.4) Co,, M ^ (2C(<i - 2, |) - (cPrM- 

The positive part of the dimension spectrum is 3}, with a = dimji^(Py), while 

the full dimension spectrum S also contains the points fl5.2p off the real line. The metric 
dimension of the spectral triple is 3. 


Proof. The spectrum of the Dirac operator on the Poincare homology sphere Ya, with the 
correct multiplicities, can be computed explicitly using the method of [2] of generating func¬ 
tions, see [5l]. In the case of the trivial spin structure, it is shown in [5l] that one can separate 
the spectrum into 60 arithmetic progressions {^ + k + 60j} with multiplicities interpolated 
by 60 explicit polynomials Pk{^ + k + 60j) = Mult(| + k + 60j), which satisfy 


59 


Y, PM 


1 2 1 
= -u — 


k=0 


This implies that 


59 


Sv.(A) = E E ^ ^ ^ 


and that 


k=0 j& 


59 


A 


Cdv, («) = ^ ^ ^ 

A:=0 jGZ 

hence (15.411 then follows as in Proposition 12.61 


l^ + k + 60j 


- 553 (A) 

120 ’ 


- Cd -i(s), 

120 sr 


□ 


Using the same technique that we used in the previous construction, we can then compute 
the leading terms in the expansion of the spectral action. In this case we have a completely 
explicit description of the poles off the real line, so we also obtain a more explicit description 
of the log oscillatory corrections to the spectral action, which in this model behaves exactly 
as in the case of Proposition 13.11 


Proposition 5.3. Let Vy be the fractal arrangement described above, to which we assign a 
spectral triple as in Proposition \5.^ The contribution to the expansion of the spectral action 
Tr(/(Dpy/A)) coming from the points of the dimension spectrum that lie on the positive real 
line is given by 


(5.5) 


(Aa) 


\Cc{Vy){^) 


120 


/s - Aa 


C/:()Py)(l) 


120 


fi + (Aa)' 


.C(^-2,|) 


K(^ 4 : 


120 log(2 + 0) 


■u 


where a = dimj 5 (Py), while the contribution to the expansion of the spectral action coming 
from the non-real points of the dimension spectrum is given by a Fourier series 


A' 


Ea 

m^O 


2, 2 ) 4C(®m! 2 ) 
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with 

POO 

fsm = / du. 

Jo 

The series converges absolutely to a smooth function of the angular variable 9 = iog° 2 +</i) "^^d 
27rZ. 

Proof. The argument is exactly as in Proposition 13.51 where we now have 2 < a < 3 and the 
residue 71^ (as well as those at the poles off the real line) is given by 

= ReSs=^CciPy){s) = iQg^2 + (f) • 

The complex poles Sm of Cc{s) he within the region of absolute convergence of the series 
dehning the Hurwitz zeta function, hence the size of the terms |C(sm ~ 2, |) — |)| is 

controlled by a term |(C((T — 2, |)| + ^|C(c’‘) |)l- The convergence of the Fourier series above 

2TTim 

is then controlled by the convergence of which converges absolutely to a 

smooth function of the periodic angle variable 6, as shown in Proposition 13.11 □ 


Corollary 5.4. Consider a (Euclidean, compactified) spacetime model of the form x Vy, 
with Vy the fractal arrangement as above, and with (3 the size of the compactification. Then 
the contribution to the expansion of the spectral action SsjxVyi-^) coming from real points 
of the dimension spectrum is given by 

^ ^ 120 ^ 120 ^ 12Olog(2 + 0) 

with a = dimuifPy) and f}i, 1)3, as in fl3.1ip and fl3.12p . The contribution to the points of 
the dimension spectrum that are off the real line is a Fourier series 


A^+i 


2/3a° 


120 log(2 + 0) 


E 27rim , j I 

(Aa)iog(2+« -) - ((sm,-))hs„„ 


m^O 


where the coefficients are given by 

f).. =2 


h{p^)p^^dp. 


The series converges to a smooth function of 6 = (Qg° 2 +^) 
Proof. The result follows exactly as in Proposition 13.61 


□ 


6. Slow roll inflation in fractal universes 

In the case of a compactihed Euclidean spacetime of the form Sj^ x Sf it was shown in 
[9] that perturbations of the Dirac operator by a scalar held i-A + cjp produce, in 
the calculation of the spectral action, a potential E((/)) for the scalar held, obtained as a 
combination of functions V((/)^/A^) and W(0^/A^) with 

POO 

IsfjJa / h{u)du 

Jo 

+7rA^l3a^V{(fyA^) + ^A^PaWif^A^), 


POO 

Tr(h((D^ + y)^)/A^))) ~ vrA'^/Sa^ / uh{u)du — — 

Jo 2 
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where the functions are of the form 



■OO 


>V(a:) = [ h{u)d- 


u{h{u + x) — h{u))du, 


\u. 


In [9] it was first observed that the potential obtained in this way has the typical shape 
of the slow-roll inflation potentials. Whether one can accommodate a satisfactory inflaton 
model within the spectral action paradigm is still debated. It was first proposed that the 
Higgs field might play the role of inflaton field, but the possibility of a Higgs based inflation 
scenario in the noncommutative geometry model was ruled out as incompatible with the 
measured value of the top quark mass in [S], based on constraints coming from the CMB 
data. In the noncommutative geometry models of gravity coupled to matter, the Higgs sector 
arises as inner fluctuation of the Dirac operator in the non-commutative hber directions of 
an almost-commutative geometry. By contrast, even in the pure gravity case, where there 
is no hnite non-commutative geometry, one can introduce a scalar perturbation of the Dirac 
operator of the kind described above, see the discussion in §1.4 of [7]. 

It was observed in [4T], |12] and [7], where the construction is generalized for spherical 
space forms and Bieberbach manifolds, that the shape of the resulting slow-roll inflation 
potential H(0) distinguishes between (almost all of) the different possible topologies and 
determines detectable signatures of the cosmic topology in the slow-roll parameters (which 
in turn determine spectral index and tensor-to-scalar ratio) and in the form of the power 
spectra for the scalar and tensor fluctuations. The result is similar in the case of the spherical 
space forms W = *S'f/r, with H(0) replaced by (#r)“^l/(0), see [7], [H]. In this section 
we discuss how the inflation potential changes in the case of x P, where V is either an 
Apollonian packing or 3-spheres or a configuration based on the fractal dodecahedron and 
the Poincare homology sphere. In particular, we show that, in models of inflation based 
on the spectral action functional, the shape of the inflation potential changes depending 
on the fractal structure, hence the potential detect measurable effects of the presence and 
the type of fractal structure. In particular, it follows that, in a model of gravity based on 
the spectral action, the presence of fractality in the spacetime structure leaves a detectable 
signature in quantities, like the slow-roll parameters and the power spectra for the scalar 
and tensor fluctuations, that are in principle measurable in the CMB, modulo the problem 
of determining the unknown parameter f3 of the model, already discussed in [TT], |42j . 


We perturb the Dirac operator D of a packing V of 3-spheres by a scalar field (p. Corre¬ 
spondingly the spectral action is modihes by terms that determine a potential for the scalar 
held. We discuss the effect of the real points of the dimension spectrum (Proposition 16.11) 
and of the huctuations coming from the points off the real line (Proposition [6^ on the shape 
of the potential. 


Proposition 6.1. Let V he either a packing of 3-spheres satisfying the three conditions of 
Definition \3.d[ or a configuration of Poincare homology 3-spheres arranged according to the 
fractal dodecahedron construction of l id, il Consider the product geometry SjjXV ofV with 
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a circle of size f3. Then the spectral action satisfies 

TiihiiD^ + fi^)/K^))) ~ 2/3 (A^ C£(3) Os - i Cz:(l) Oi) 

+2/3A-+1 (C(a - 2,1) - iC(c^, §)) f). 

(6.1) +vrA4/3C£(3) W{ctfi/K^) 

+4/3A‘^+i (C(cr - 2 , 1) - lC(c^, D) 

+‘5r^(A), 

where the last term collects the fluctuations coming from the log-oscillatory terms contributed 
by the poles of the zeta functions that are off the real line, described in Proposition ^.^ below. 
The potentials V, W, and Ua are, respectively, given by 


( 6 . 2 ) 


V(x) = / u{h{ux) — h{u))du, W{x) = / h{u)du, 


(6.3) 


Ua{x) = / x) — h{u))du, 


where a = crfiV) is the packing constant fl2.5l) . TZa is the residue of (c{s) at s = a, and Oi, 
Os, Oa are as in 03.111) . 03.121) . 

Proof. The argument follows directly from Proposition 13.61 along the lines of Theorem 7 of 
[9]. We have 

roo 2 

/ h{p^)p^dp = 77 / uh{u)du, 

Jo 2 Jq 

which gives rise to the term 

noo 

/ u{h{ux) — h{u))du 

Jo 

in the V part of the potential; similarly, we have 

noo 2 

/ h{p^)pdp = - h{u)du, 


which gives the term 


h{u)du 


in the W part of the potential, and the term 


^OO 1 


which gives the term 


in the U„ part of the potential. 


h{p )p'^dp = - u 2 h{u)du, 
^ .In 


u 2 (h{ux) — h{u))du 


□ 
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Proposition 6.2. Under the assumptions of Definition \3.^ the fluctuation terms 
are of the form (A) = where Ssiy^-p{A)°‘^'^ is as in fl3.13p . Let 

0 osc( 0 )^^ 6 e the approximation to that only counts the contribution of non-real poles 

with |9'(s)| < R. The oscillatory term is approximated by a sequence 

(6.4) 


4/?^ A” 


7+1 




(<^r3 


log bn 


’(C(* 


n,j,m 2 , 


1 

4 


C{^n,j,my 



Res 


where n —)■ cxo and R oo, with 

/ OO 

+ x) - h{u))du, 

and where Sn,j,m = (^n,j + o,re the complex zeros of the series of self-similar 

strings Cn with the lattice property approximating the complex poles of(c{P){s). 


Proof. The result is obtained as in Proposition 16.11 using the results of Proposition I3.6[ □ 


The slow-roll inflation potential V{4>) is obtained from the combination of functions V, 
W, U that appears in the expansion of the spectral action above. 

Corollary 6.3. The function 

TiA^flCdS) V(</>VA^) + ^A^/^Cai) >V(0VA^) 

+4/?A-’+‘ C.T - 2, |) - |)) 'RM4>gP) + WTiZ) 

depends explicitly on the presence of fractality, through the coefficients C/:( 3)7 C/:(l)? fhe 
residue TZ^, the packing constant a, and the oscillatory fluctuations. 


Corollary 6.4. In the case of the fractal arrangement Vy of dodecahedral spaces considered 
in the previous section, the form of the fluctuations in the inflation potential is simpler and 
given by the Fourier series 


ur(0) 


^ / 27rim 3 

log(2 + 0) ^ V ^ ^ log(2 + 0) ’2 


1 2Trim 3\ 0^ 

4^^"^ ^ log(2 + 0) ’ 27 


with 


U. 


(X = 


^t(°'+iog( 2 + 0 ) -\- x) — h{u)) du. 


Proof. This follows directly from the results of the previous section, by proceeding as in the 
proof of Proposition 16.11 □ 


7. Conclusions and further questions 

7.1. Conclusions. In this paper we considered a model of gravity based on the spectral ac¬ 
tion functional. This is known to recover, via its asymptotic expansion, the usual Einstein- 
Hilbert action with cosmological term, along with modihed gravity terms (conformal and 
Gauss-Bonnet gravity). We considered simple models of (Euclidean, compactihed) space- 
times of the form x V, where fl is the size of the R^-compactihcation and P is a fractal 
conhguration built out of 3-spheres (Apollonian packings) or of other spherical space forms 
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(Sierpinski fractals). We evaluated the leading terms of the expansion of the spectral action, 
using information on the zeta function of the Dirac operator of a spectral triple, and we 
compared them, respectively, with the corresponding terms in the simpler case of x S'f 
(spatial sections given by a single sphere of radius a) or xY where F is a spherical space 
form, in particular the Poincare homology sphere (dodecahedral space). We regard the case 
of Sjj X V, where V is an Apollonian packing of 3-spheres or a conhguration obtained from 
a Sierpinski fractal dodecahedron, as a model of possible presence of fractality in spacetime 
geometry: a version of Packed Swiss Cheese Cosmology models. We showed that the result¬ 
ing leading terms of the expansion of the spectral action for xV differ from those of the 
ordinary x (or x Y) case in the following ways: 

• The term 2A'^/3a^()3 —respectively corresponding to the cosmological and the 
Einstein-Hilbert term, are replaced by terms of the form 2A'^/3C£(3)[)3 — (3C,ciX)^i.i 
which can be seen as a zeta regularization of the divergent series of the 3-sphere terms 
summed over the packing. 

• There is an additional term in the gravity action functional of the form 

(^C(<^ - 2, |) - |)) UA, 

where C('5,2:) is the Hurwitz zeta function and a is the packing constant of P, with 
Po- the residue at a of the zeta function Cc{s) of the fractal-packing. For a fractal 
dodecahedron a is the Hausdorff dimension, while for an Apollonian packing it is 
conjecturally the Hausdorff dimension of the residual set. In both cases this term 
detects modihcations to the gravity action functional due to the presence of a fractal 
structure. 

• This additional term is further corrected by a Fourier series of log-oscillatory fluc¬ 
tuations, coming from the presence of points of the dimension spectrum off the real 
line (another purely fractal phenomenon). 

• The perturbation i-A -\- (jY of the Dirac operator determines a slow-roll inflation 
potential V (0) for the held 0. The shape of the potential detects the presence of 
fractality through the coefficients C/:(3)) C/:(l)) fhe packing constant a, and the residue 
at a of C/:('S), and oscillatory huctuations. 

7.2. Further questions. There are a number of questions, both mathematical and physical, 
that arise in relation to improving the model of gravity, based on the spectral action, on 
cosmologies exhibiting fractality. On the mathematical side, as we have seen, one needs a 
better understanding of the properties of higher-dimensional Apollonian packings, especially 
with respect to characterizing the presence of exact self-similarity, extending results like [T5] 
beyond the 2-dimensional case of circle packings. 

From the physics viewpoint, the logic we followed in this paper is along the lines of sev¬ 
eral other recent results, where one considers a classical model of (Euclidean, compactified) 
spacetime and computes what the expansion of the spectral action looks like, either with the 
full infinite series of the asymptotic expansion, or at least with the leading terms up to order 
A°, see for instance [9], [10], [25], [27|, [5l]. Under this perspective, one can consider the 
classical Packed Swiss Cheese Cosmologies. These originate in a spacetime model introduced 
in the late ’60s in m as a prototype model of isotropic but non-homogeneous spacetimes 
(as opposed to, for instance, the Bianchi IX examples of homogeneous non-isotropic spaces). 
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In other words, the original constrnction of the Packed Swiss Cheese Cosmology is dictated 
by imposing certain kind of regnlarity reqnirements on the geometry. These classical models 
have then been tested as possible models of fractal and mnltifractal strnctnres in spacetime. 
Here we investigate how the spectral action fnnctional, seen as our choice of action functional 
of gravity, behaves on this specihc classical geometry. A hrst important question in this di¬ 
rection, which would make the model more realistic (closer to the original construction of 
im) would be to start from a (Euclidean, compactihed) Robertson-Walker spacetime and 
carve out balls, so as to obtain a residual Apollonian sphere packing, rather than adopting 
the simplihed model we considered here of a product V x S^. In terms of computations of 
the spectral action, this would mean adapting the computation for the Roberston-Walker 
metric of [10] to the resulting fractal packing, rather than (as we did here) adapting the 
computations of [9] for to the case of the Apollonian packings of 3-spheres, or the fractal 
dodecahedral packing of Poincare homology spheres. A second question would be to adopt 
a different viewpoint and derive spacetime models (possibly with some form of fractal struc¬ 
tures or of noncommutativity) from a least action principle applied to the spectral action 
functional. More explicitly, the question would be whether such a variational principle can 
be reinterpreted in terms of classical theory as a modihed gravity model, with an effective 
stress-energy tensor (as for instance in the case of /(i?)-modihed gravity). 
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